REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis 
Highway,  Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a 
collection  of  information  if  it  does  not  display  a  currently  valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 

1 .  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE 

12-06-2007  Journal  Article 

3.  DATES  COVERED  (From  -  To) 

4.  TITLE  AND  SUBTITLE 

The  Impact  of  Manifold-to-Orifice  Turning  Angle  on  Sharp-Edge  Orifice  Flow 
Characteristics  in  both  Cavitation  and  Non-Cavitation  Turbulent  Flow  Regimes 
(Preprint) 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

W.H.  Nurick  &  T.  Ohanian  (Science  and  Technology  Applications  LLC); 

D.G.  Talley  &  P.A.  Strakey  (AFRL/PRSA) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

50260538 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Air  Force  Research  Laboratory  (AFMC) 

AFRL/PRSA 

10  E.  Saturn  Blvd. 

Edwards  AFB  CA  93524-7680 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

AFRL-PR-ED-JA-2007-325 

9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR’S 
ACRONYM(S) 

Air  Force  Research  Laboratory  (AFMC) 
AFRL/PRS 
5  Pollux  Drive 

Edwards  AFB  CA  93524-7048 


11.  SPONSOR/MONITOR’S 
NUMBER(S) 

AFRL-PR-ED-JA-2007-325 


12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 


Approved  for  public  release;  distribution  unlimited  (PA  #07234 A). 


13.  SUPPLEMENTARY  NOTES 

Submitted  for  publication  in  Journal  of  Spacecraft  and  Rockets.. 

14.  ABSTRACT 

The  approach  taken  was  to  analyze  the  results  in  a  manner  consistent  with  application  by  design  engineers  to  new  and  existing 
applications,  while  providing  some  insight  into  the  processes  that  are  occurring.  This  paper  deals  with  predicting  the  initiation  of  cavitation, 
cavitation  impacts  on  Cc,  and  non-cavitation  impacts  on  Cd  from  L/D  of  5  sharp-edge  orifices  with  both  single  angle  and  compound  angle 
directional  flow.  The  results  show  that  in  the  cavitation  regime,  Cc  is  controlled  by  the  cavitation  number,  where  the  data  follows  the  power 
with  Kcav,  and  inception  of  cavitation  occurs  at  a  Kcav  of  1.8.  In  the  non-cavitation  regime  Cd  is  controlled  by  Reynolds  number,  and  the  head 
loss  coefficient,  KL  for  all  angles  is  a  function  of  the  manifold-to-orifice  velocity  ratio.  Compound  angle  orifices  Cc  and  KL  were  found  to  be 
influenced  more  by  the  initial  turning  angle  than  the  orifice  turning  angle.  In  the  non-cavitation  regime  for  conditions  where  the  cross  velocity  is 
0  the  data  are  consistent  with  the  first  order  equation  relating  HL  to  the  dynamic  pressure  where  KL  is  constant,  which  is  consistent  with  in-line 
orifices. 


15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

18.  NUMBER 

19a.  NAME  OF  RESPONSIBLE 

OF  ABSTRACT 

OF  PAGES 

PERSON 

Dr.  Douglas  Talley 

a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

19b.  TELEPHONE  NUMBER 

SAR 

21 

(include  area  code) 

Unclassified 

Unclassified 

Unclassified 

N/A 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  239.18 


THE  IMPACT  OF  MANIFOLD-TO-ORIFICE  TURNING  ANGLE  ON  SHARP- 
EDGE  ORIFICE  FLOW  CHARACTERISTICS  IN  BOTH  CAVITATION  AND 
NON-CAVITATION  TURBULENT  FLOW  REGIMES  (preprint) 


By: 

W.  H.  Nurick  &  T.  Ohanian 
Science  &  Technology  Applications  LLC  (STA) 
Moorpark,  CA 

And 


D.G.  Talley  &  P.  A.  Strakey 
Air  Force  Research  Laboratory 
Edwards  Air  Force  Base  CA 


ABSTRACT 

The  approach  taken  was  to  analyze  the  results  in  a  manner  consistent  with  application  by 
design  engineers  to  new  and  existing  applications,  while  providing  some  insight  into  the 
processes  that  are  occurring.  This  paper  deals  with  predicting  the  initiation  of  cavitation, 
cavitation  impacts  on  Cc,  and  non-cavitation  impacts  on  Cd  from  L/D  of  5  sharp-edge 
orifices  with  both  single  angle  and  compound  angle  directional  flow.  The  results  show 
that  in  the  cavitation  regime,  Cc  is  controlled  by  the  cavitation  number,  where  the  data 
follows  the  V2  power  with  Kcav,  and  inception  of  cavitation  occurs  at  a  Kcav  of  1.8.  In 
the  non-cavitation  regime  Cd  is  controlled  by  Reynolds  number,  and  the  head  loss 
coefficient,  KL  for  all  angles  is  a  function  of  the  manifold-to-orifice  velocity  ratio. 
Compound  angle  orifices  Cc  and  KL  were  found  to  be  influenced  more  by  the  initial 
turning  angle  than  the  orifice  turning  angle.  In  the  non-cavitation  regime  for  conditions 
where  the  cross  velocity  is  0  the  data  are  consistent  with  the  first  order  equation  relating 
Hl  to  the  dynamic  pressure  where  KL  is  constant,  which  is  consistent  with  in-line  orifices. 
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INTRODUCTION 


For  application  to  rocket  engine  injectors,  cavitation  can  occur  in  the  entrance  from  flow 
passages  into  sharp-edge  orifices.  The  actual  onset  of  cavitation  depends  on  the  orifice 
edge  sharpness,  L/D,  upstream  pressure,  cross  velocity,  orifice  entrance  angle,  and  back 
pressure.  In  application  where  the  engine  is  “deep”  throttled  cavitation  can  occur  if  the 
chamber  pressure  is  sufficiently  low  and/or  the  upstream  flow  causes  a 
contraction/expansion  to  support  separation.  The  inability  to  predict  the  occurrence  of 
cavitation  can  lead  to  a  catastrophic  failure  of  an  engine  and  loss  of  mission. 

In  1976  Nurick1  published  a  paper  discussing  cavitation  for  in-line  sharp-edge  orifices 
and  proposed  a  simple  1st  order  model  for  predicting  cavitation.  Photographic  evidence 
was  also  presented  showing  that  after  the  inception  of  cavitation,  the  vena-contracta 
reattachment  point  started  to  move  downstream  toward  the  exit  of  the  orifice.  This  also 
resulted  in  a  decrease  in  discharge  coefficient  (Cd)  until  it  reached  the  separation  value.  In 
some  cases  the  orifice  remained  attached  even  at  a  cavitation  number  of  1  while  in  other 
cases  hydraulic  flip  occurred  in  the  region  of  Kcav  of  1.3  or  less.  In  2007  Nurick, 
Ohanian,  Talley  and  Strakey2  (to  be  published)  wrote  a  paper  specifically  dealing  with 
sharp  edged  orifices  having  a  90  degree  orientation  relative  to  the  manifold  feed.  This 
study  revealed  that  the  simple  relationship  could  be  extended  to  90  degree  orifices  in 
cross  flow  although  the  Cc  and  KL  relationships  needed  to  account  for  the  manifold  to 
orifice  velocity  ratio. 

The  early  Nurick1  results  were  subsequently  utilized  in  other  studies3’  4’ 5  that  included 
CFD  modeling  to  extend  model  applicability  and  in  some  cases  verify  codes.  In  the  last 
10-15  years  CFD  modeling  has  improved  considerably  and  predictions  are  now  being 
made  that  show  the  impact  of  key  design/operating  variables  on  both  the  discharge 
coefficient  (Cd)  and  contraction  coefficient  (Cc).  Unfortunately,  there  is  little  if  any  test 
data  that  can  be  utilized  to  verify  these  predictions.  It  is  the  hope  of  the  authors  that  the 
results  presented  in  this  paper  will  be  useful  in  that  goal. 

This  paper  is  an  extension  of  the  earlier  work1,  2  and  includes  varying  the  turning  angle 
from  as  low  as  60°  up  to  120°,  for  an  orifice  L/D  of  5.  The  approach  taken  in  this  paper  is 
to  present  the  results  in  a  manner  consistent  with  application  by  design  engineers  to  new 
and  existing  applications,  while  providing  some  insight  into  the  processes  that  are 
occurring.  This  paper  deals  with  predicting  the  initiation  of  cavitation  and  the  impact  of 
the  variables  listed  above  on  Cc,  Cd,  and  KL.  In  addition  the  impact  of  the  same  variables 
in  the  turbulent  flow  regime  is  also  discussed  and  correlations  are  provided. 

TEST  FACILITY  AND  TEST  SETUP 


FACILITY  DESIGN 

The  experimental  investigation  was  carried  out  at  the  Air  Force  Research  Laboratory 
cold-flow  injector  characterization  facility,  a  simplified  schematic  of  which  is  shown  in 
Figure  1.  Water,  which  is  used  as  a  simulant  for  liquid  oxygen,  is  stored  and  pressurized 
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in  a  35  ft3  tank.  The  injector  inlet  flow  rate  is  controlled  with  a  throttling  valve  and 
measured  with  one  of  several  turbine  flow  meters  arranged  in  parallel  to  cover  a  wide 
range  of  flow  rates.  Downstream  of  the  injector,  another  turbine  flow  meter  measures  the 
outlet  flow  rate  and  a  back-pressure  regulator  maintains  the  injector  fluid  pressure.  Fluid 
pressures  of  up  to  1500  lb/in  and  manifold  velocities  up  to  130  f/s  may  be  produced  in 
this  manner.  The  injector,  as  illustrated  in  Figure  2,  consists  of  an  interchangeable 
stainless  steel  plate  with  a  sharp-edged  hole  machined  in  the  center  of  the  plate.  This  is 
mounted  against  the  fluid  manifold,  which  contains  a  0.25  in  square  channel  machined 
into  the  manifold.  A  plenum  and  screen  at  the  inlet  and  exit  of  the  manifold  reduces  the 
velocity  before  the  flow  enters  the  manifold  in  order  to  generate  a  reproducible  flow  field 
at  the  entrance  to  the  orifice.  The  distance  from  the  manifold  inlet  screen  to  the  orifice 
entrance  is  2.5  in  or  9  manifold  channel  widths.  The  entire  arrangement  is  secured  inside 
an  optically  accessible  pressure  chamber  which  is  rated  to  2000  psi. 


Pump 

Figure  1  -  Schematic  of  AFRL  Cold  Flow  Test  Facility 

The  orifices  are  either  pilot  drilled  and  reamed,  or  made  by  electrical  discharge 
machining  with  a  diameter  tolerance  of  +/-  13  pm  and  inlet  edge  radius  to  orifice 
diameter  ratio  of  less  than  .003,  ensuring  a  sharp  edged  inlet.  Chamber  pressure,  orifice 
pressure-drop,  and  inlet  and  outlet  flow  rates  are  recorded  by  a  12  bit  analog  to  digital 
conversion  board  and  the  data  is  stored  on  a  personal  computer.  Experiments  are 
typically  conducted  by  setting  the  fluid  pressure  and  flow  rates  to  a  predetermined  value, 
with  the  chamber  pressure  being  gradually  increased  while  the  data  acquisition  system 
records  flow  rates  and  pressures.  This  allows  for  a  large  amount  of  data  to  be  collected  in 
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a  relatively  short  period  of  time.  Orifice  pressure  drop  and  chamber  gas  pressure  are 
measured  within  an  accuracy  of  +/-  0.25%.  Manifold  velocities  are  held  constant  to 
within  +/- 1.5%  during  the  experiment.  Because  of  the  difficulty  associated  with  directly 

Transducer 

Port 


Figure  2  -  Injector  Schematic 

measuring  the  orifice  flow  rate  inside  of  the  pressurized  vessel,  the  orifice  flow  rate  is 
measured  by  subtracting  the  manifold  outlet  flow  rate  from  the  manifold  inlet  flow  rate. 
The  experimental  error  associated  with  the  discharge  coefficient  measurement  is  limited 
by  the  accuracy  of  the  flow  meters  which  is  +/-  0.5%.  This  translates  to  an  error  on  the 
discharge  coefficient  data  of  about  +/-  0.5%  at  the  lowest  manifold  flow  rates  and  highest 
orifice  flow  rates,  to  +/-  10%  at  the  lowest  orifice  flow  rates  and  highest  manifold  flows. 
A  typical  error  for  the  intermediate  flow  rates  is  on  the  order  of  +/-  4%. 

MANIFOLD  CONFIGURATION 

For  the  data  in  this  paper  the  manifold  was  operated  such  that  fluid  enters  the  inlet  port 
and  a  portion  of  the  flow  exits  the  manifold  and  the  remainder  enters  the  orifice.  A 
sketch  of  the  configuration  is  shown  in  Figure  3. 
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TEST  CONFIGURATION 


Flow  In 


Orifice  Flow  Out 


Figure  3  -  Manifold  /  Orifice  Entrance  Test  Configuration 
ORIFICE  CONFIGURATIONS 

Both  single  and  compound  angle  orifice  designs  were  tested.  As  shown  in  Figure  4 
below,  for  the  single  angle  direction  feed,  the  orifice  is  oriented  such  that  the  flow  enters 
the  orifice  by  turning  in  the  direction  of  its  angle;  while  for  the  compound  angle  direction 
flow  the  flow  must  first  turn  90°  then  immediately  turn  in  the  direction  of  the  orifice 
angle. 


Single  Angle  Direction  Feed 


W2J  V2,P2 


Compound  Angle  Direction  Feed 

v2 

A 


Figure  4  -  Illustration  of  Flow  Turning  For  Both  Single  &  Compound  Angle  Orifices 

The  range  of  operating  conditions  and  orifice  geometries  typical  of  liquid  rocket  injectors 
was  studied  and  are  given  in  Table  1. 
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Table  1  -  Test  Variables  and  Range  of  Testing 


Orifice  Diameter  (in) 

Length  to  Diameter  Ratio 
Fluid  Pressure  (psi) 

Back-Pressure  (psi) 

Cross- Velocity  (f/s) 

Orifice  Red 
Manifold  Rew 

Manifold  Dynamic  Head  to  Orifice  Static  Dp  Ratio 


0.0468  -  .0839 
3  -  10 
100  -  1500 
13.05  -  1500 
0  -  60.7 

l.OxlO4-  3.0xl05 
6.0xl03 

5.0xl0'5- 6.0x10° 


ANALYSIS  APPROACH 


CAVITATION  REGIME  (TURBULENT  FLOW) 


The  definition  of  cavitation  number  should  take  into  account  all  the  dynamic  variables 
that  impact  the  process.  For  high  Re  number  flow  Cc  is  also  a  function  of  A2/Ai  and 
therefore  the  cavitation  number  should  include  the  manifold  contraction,  any  flow  turning 
as  well  as  the  processes  to  the  exit.  Consistent  with  this  the  definition  of  cavitation 
number  utilized  in  this  study  is  stated  as: 


K 


cav 


$  ~  Py 

Pl~P2 


(1) 


The  derivation  of  relationships  between  Cd  and  Cc  depend  on  various  assumptions  unique 
to  the  flow  characteristics  of  the  configuration  tested.  For  the  tests  in  this  program  the 
manifold  velocity  can  not  be  neglected  in  the  derivation  of  discharge  coefficient.  For  this 
case  the  discharge  coefficient  is  defined  as: 

cj  ■  r - z — 1  ; — 7T7  <2> 

Combining  Equations  1  and  2  in  addition  to  the  Bernoulli  equation  with  the  assumption 
that  the  head  loss  between  the  manifold  and  the  vena-contracta  is  small  results  in 
Equation  3: 


CK. 


1/2 


(VA 

v22  1 

UJ 

W 

-.1/2 


c.2(£„,-i) 


(3) 


For  conditions  where  all  the  manifold  flow  enters  the  orifice  Equation  3  reduces  to: 
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Q  = 


CK, 


1/2 


-1/2 


i+f^  |  cAk„,-i) 

\  A1 


And  for  in-line  orifices  where  Vi  «  V2  and  Vout  =  0: 

cAkJ' 


1/2 


Q  = 


-.1/2 


Cc2 


(4) 


(5) 


Then  when  the  product  of  area  ratio  times  Cc  squared  is  small  with  respect  to  1  this 
reduces  to: 


Cd  =  CcKj'2  (6) 

This  expression  was  used  in  the  early  Nurick1  study  and  subsequently  adopted  by  others. 

NON-CAVITATION  REGIME  (TURBULENT  FLOW) 

In  basic  fluid  dynamics  texts,  flow  around  bends  is  defined  by  a  loss  coefficient  which  is 
included  in  the  specification  of  overall  head  loss.  Consistent  with  this  tradition  the  data 
in  the  non-cavitation  turbulent  flow  is  both  analyzed  and  correlated  in  terms  of  the  impact 
on  overall  head  loss.  Based  on  constant  velocity  studies  where  turning  angle  has  been 
varied  the  loss  coefficient  has  been  found  to  be  a  function  of  the  turning  angle. 

Therefore,  the  loss  coefficient  (KL)  is  defined  by: 

Hl  -  Kl(6)£-V?  (7) 

2  g 


For  conditions  where  the  velocity  is  not  constant  then  Kl  should  also  be  a  function  of 
V1/V2  or: 


hl  =  kl 


ev^ 


v. 


—V22 

2  g 


(8) 


The  Bernoulli  equation  was  used  to  define  HL: 


HL=(Pl-P2)  +  ^-(Vl2-V22)  (9) 

2g 
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RESULTS  AND  DISCUSSION 


CAVITATION  REGIME 

It  is  important  to  first  define  the  flow  regimes  so  that  the  data  can  be  related  to  specific 
processes.  There  are  three  flow  regimes  of  interest  in  the  cavitation  regime.  The  first  is 
the  onset  of  cavitation  (inception  of  cavitation)  where  bubbles  are  formed  at  the  orifice 
entrance  and  continue  to  flow  downstream  until  they  are  converted  back  to  liquid  in  the 
recovery  area.  The  second  is  full  cavitation  where  the  bubbles  coalesce  forming  a  vapor 
cavity  whose  length  varies  as  the  cavitation  number  is  lowered.  The  third  is 
supercavitation  where  the  vapor  pocket  attachment  has  moved  to  the  orifice  exit  and 
beyond  but  the  flow  still  acts  as  attached.  For  this  study  only  the  90  degree  orifice  angle 
configuration  experienced  supercavitation.  For  all  other  angles  it  did  not  appear  to  occur. 
The  physical  reason  is  unknown  at  this  time  and  its  determination  would  require 
additional  measurements  including  photographic  evidence  of  the  vapor  cavity 
characteristics. 

All  three  of  these  regimes,  in  addition  to  separation,  are  illustrated  in  Figure  5.  Note  that 
hydraulic  flip  can  occur  at  any  point  in  the  cavitation  regime  depending  on  the  orifice 
L/D  as  well  as  flow  conditions.  There  also  appears  to  be  potential  for  “hysteresis”  to  also 
occur  at  start-up  as  indicated  by  the  two  identical  run  conditions  shown  in  Figure  5.  Note 
that  in  one  case  the  orifice  was  separated  at  the  outset  but  at  about  the  point  were 
cavitation  inception  occurs  it  immediately  flips  to  the  non-cavitation  condition.  In  the 
other  test  it  appears  that  the  flow  is  attached  at  the  onset  but  not  in  separation  since  the 
data  curves  upward  until  about  a  Kcav  of  1.17  where  Cc  becomes  constant  and  then,  as  is 
normal,  transitions  to  the  non-cavitation  condition.  The  transition  zone  is  however 
impacted  by  the  ability  of  the  flow  to  remain  attached  and  other  processes  most  likely 
occurring  at  the  orifice  exit.  These  resulting  impacts  on  the  cavitation  characteristics  are 
noted  in  Figure  5  by  variable  Cc.  Also  shown  is  the  constant  Cc  condition  here  the  flow 
smoothly  transitions  from  non-cavitation  to  cavitation  and  remains  attached.  This  later 
condition  is  the  basis  for  the  correlation  presented  in  this  paper 


Figure  5  -  Sharp-Edge  Orifice  Data  Depicting  All  Flow  Regimes 
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Single  Angle  Direction 

Inception  of  Cavitation  -  The  inception  of  cavitation  occurs  where  the  discharge 
coefficient  either  abruptly  changes  to  a  decreasing  value  with  decreasing  cavitation 
number  or  more  smoothly  transitions,  both  resulting  in  a  constant  slope  consistent  with 
Equation  3.  Both  characteristics  are  illustrated  in  Figure  6.  Interestingly  the  cavitation 
data  all  follow  the  relationship  of  Equation  3  suggesting  that  the  denominator  for  a  given 
test  is  constant.  This  is  not  unexpected  since  the  tests  were  run  such  that  the  velocity 
ratio  (V1/V2)  remained  constant  for  a  given  test. 
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Figure  6  -  Cavitation  Number  Vs  Discharge  Coefficient 

For  each  test  sequence,  where  the  upstream  conditions  were  set  and  the  back  pressure 
systematically  lowered,  the  data  were  analyzed  to  determine  the  point  where  the  data 
deviated  from  the  straight  line  relationship.  This  is  also  the  point  where  the  flow  is  no 
longer  choked  and  transitions  into  the  non-cavitation  turbulent  regime.  Analysis  of  the 
data  revealed  that  the  cavitation  number  where  cavitation  inception  begins  remains 
constant  (within  experimental  error)  at  all  conditions  regardless  of  turning  angle  or 
upstream  conditions  at  Kcav  ~  1.8.  In  some  cases  there  is  a  transition  region  where  Cc 
varies  until  the  flow  fully  chokes  and  Cc  remains  constant.  When  this  occurs  Cc  becomes 
constant  at  Kcav  between  ~  1.5- 1.6. 

Full  Cavitation  -  Since  the  impact  of  upstream  conditions  on  the  contraction  coefficient 
should  be  related  to  the  non-symmetric  nature  of  the  flow  due  to  acceleration  from 
bending  as  well  as  both  area  reduction  and  attachment  on  the  far  wall.  These 
characteristics  should  impact  the  vena-contracta  formation  process.  Equation  3  can  be 
solved  for  Cc  and  the  result  is: 
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Therefore  it  is  expected  that  Cc  will  be  a  function  of  the  velocity  ratio  (area  ratio  is 
constant  for  these  tests  and  the  ratio  Vout/V2  is  dependent  on  V1/V2).  For  different  turning 
angles  it  would  be  expected  that  Cc  would  also  be  a  function  of  angle.  The  contraction 
coefficient  is  therefore  a  function  of: 


Ce  =  (13) 

For  a  given  orifice  the  tests  were  conducted  at  a  constant  velocity  ratio,  and  since  the  area 
ratio  is  small  it  would  be  expected  that  the  contraction  ratio  would  follow  the  Equation  6 
relationship.  To  test  this  hypothesis  all  data  were  first  evaluated  according  to  Equation  6 
and  typical  results  are  presented  in  Figure  7. 


k„.,.4-s 

Figure  7  -  Typical  Test  Results  Plotted  According  to  Equation  6 


These  results  and  in  fact  all  test  data  show  the  same  trend  that  confirm  that  the 
contraction  ratio  is  a  function  of  both  turning  angle  and  velocity  ratio. 
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Even  though  it  would  be  expected  that  using  Equation  4  or  Equation  6  would  yield  near 
identical  results  the  test  data  was  reduced  to  define  Cc  using  Equation  10  for  each  turning 
angle.  The  results  are  plotted  in  Figure  8  as  a  function  of  V 1/V9.  For  the  60°  turning 
angle  the  results  show  that  Cc  initially  increases  with  velocity  ratio  then  decreases  similar 
to  that  of  the  other  angles.  As  the  manifold  exit  flow  rate  approaches  zero  the  Cc  appears 
to  “sense”  the  resistance  of  a  wall  being  formed  and  therefore  seek  the  value  for  no  exit 
flow. 


Figure  8  -  Impact  of  Cross  Velocity  and  Turning  Angle  on  Cc 


The  significant  impact  of  turning  angle  and  velocity  ratio  on  Cc  suggests  that  increasing 
either  variable  results  in  a  decrease  in  the  vena-contracta  area.  This  further  suggests  that 
the  ability  of  the  accelerating  flow  to  reach  vapor  pressure  increases  as  the  flow  turning 
angle  and/or  the  velocity  ratio  increases.  This  decrease  in  the  vena-contracta  area 
increases  the  vapor  area.  A  hypothesis  for  this  phenomenon  is  that  it  is  related  to  the 
dynamics  of  the  flow  impacting  on  the  far  wall  creating  a  vapor  pocket  only  on  the  near 
wall.  Best-fit  curves  were  fit  to  each  angle  in  the  form: 


C  =  A 


(v  V  v 

1  '  +B^  +  C 


\VU 


(14) 


The  constants  for  each  angle  are  given  in  Table  2 
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Table  2  -  Constants  for  Equations  14 


Angle 

A 

B 

C 

R2 

60 

-0.5418 

0.3829 

0.5931 

0.8717 

75 

-0.1508 

-0.0180 

0.6403 

0.9458 

90 

-0.1846 

0.0202 

0.5940 

0.9211 

105 

0.0155 

-0.2506 

0.6447 

0.9945 

120 

0.0825 

-0.3034 

0.6250 

0.9916 

Application  of  this  equation  is  only  valid  within  the  limits  of  the  data  as  indicated  in 
Figure  8. 

Utilizing  the  results  from  the  cavitation  study  for  analysis  and  design  is  straight  forward. 
For  example,  given  an  orifice  and  manifold  design  as  well  as  the  expected  operating 
conditions,  Kcav  can  be  calculated  using  Equation  1.  The  value  of  Kcav  can  then  be 
compared  with  that  required  to  achieve  incipient  cavitation  (i.e.  1.8)  and  if  equal  to  or 
less  than  this  value  operation  will  be  in  the  cavitation  regime.  Then  based  on  the  design 
variables  such  as;  (1)  manifold/orifice  area  ratio,  (2)  flow  rate  ratio  and  (3)  turning  angle, 
Cc  can  be  calculated  using  Equation  14.  The  orifice  discharge  coefficient  can  then  be 
calculated  from  Equation  10. 

It  is  important  to  note  that  since  the  area  ratio  was  not  varied,  the  equations  have  only 
been  validated  at  the  ratio  A2/Ai  used  in  this  study.  However,  the  area  ratio  of  this  study 
was  sufficiently  small  that  it  could  be  neglected.  Therefore  Equation  14  has  only  been 
validated  when  the  area  ratio  can  be  neglected.  For  this  case  Equation  6  could  also  be 
used  to  define  Cd. 

Compound  Angle  Direction  Feed 

Inception  of  Cavitation  -  The  designation  of  90_xx  for  the  compound  angle  orifices 
designates  the  flow  must  first  turn  90  degrees  and  then  turn  in  the  direction  of  the  orifice 
angle  (xx). 

The  compound  angle  direction  feed  data,  Figure  9,  followed  the  identical  trends  as  the 
single  angle  direction  feed  data.  Within  the  experimental  accuracy,  cavitation  inception 
was  between  1.7  and  1.8  as  shown  in  Figure  9  (Kcav  =  1.3  to  1.34).  The  inception  of 
cavitation  for  both  the  90_60  and  90_75  degree  data  again  is  taken  at  the  beginning  of  the 
rollover  from  the  turbulent  regime  to  the  cavitation  regime.  It  should  be  noted  that  full 
cavitation  is  not  achieved  until  ~  Kcav  =1.3.  The  sharp  decline  in  Cd  in  the  non¬ 
cavitation  regime  is  due  to  the  definition  of  Cd  that  includes  V). 
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Figure  9  -  Cd  Vs  Cavitation  Number  for  Compound  Angle  90_60  Deg. 

Full  Cavitation  -  The  variation  in  the  contraction  coefficient  with  velocity  ratio  for  the 
compound  angle  direction  feed  is  shown  in  Figure  10.  The  data  are  compared  with  the 
single  angle  direction  feed  data  of  60  and  75  degrees.  Note  that  both  the  90_60  and  the 
90_75  data  nearly  fit  to  a  single  curve.  In  addition  note  that  as  the  manifold  exit  flow  rate 
approaches  zero  the  Cc  tends  to  approach  a  constant  at  a  value  equal  to  that  determined 
for  the  approach  velocity  case. 

Since  there  are  only  two  compound  angle  configurations  tested  individual  equations  are 
provided.  The  equation  is: 


C  =A 


(V  \  V 
1  >  +B^  +  C 


\V>, 


V, 


(15) 
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The  constants  for  each  configuration  are  provided  in  Table  3. 

Table  3  -  Constants  for  Equation  15 


Angle 

A 

B 

c 

90_60 

-0.2125 

-0.0.114 

0.6324 

90_75 

-0.1058 

-0.0771 

0.6259 

NON-CAVITATION  REGIME  (ATTACHED  TURBULENT  FLOW) 

The  results  for  the  non-cavitation  turbulent  flow  regime  are  divided  into  2  areas.  The 
first  section  deals  with  single  angle  direction  feed  and  the  last  section  compound  angle 
orifice  design. 

Single  Angle  Direction 

Feed  Impact  of  Orifice  Variables  on  Kl-  Since  for  each  test  the  velocity  ratio  was  held 
constant  the  data  were  initially  analyzed  using  Equation  7.  In  all  cases  the  results  were 
consistent  with  that  shown  in  Figure  11.  For  convenience  pV22/2g  in  Equation  7  is 
termed  the  dynamic  head. 


14 


400 


Dynamic  Head  -  psi 

Figure  11  -  Illustration  of  the  Linearity  of  Kl  with  the  Dynamic  Head 

For  each  test  KL  was  determined  using  Equation  7  utilizing  the  experimental  data  and  the 
results  are  shown  in  Figure  12.  The  data  fall  in  distinct  patterns  with  KL  increasing  with 
turning  angle  as  well  as  manifold  to  orifice  velocity  ratio.  This  is  consistent  with 
expectations  in  that  as  (1)  the  manifold  velocity  increases  the  acceleration  forces  increase 
then  energy  losses  should  increase  and  (2)  increasing  turning  angle  should  increase  the 
turning  loss  coefficient,  contraction  coefficient,  as  well  as  the  expansion  losses.  Since  the 
data  for  each  angle  does  not  fit  an  exponential  or  log  function  the  best  fit  equations  for 
each  angle  tested  are  provided  in  Table  4  for  convenience  of  analysis. 

Of  particular  interest  is  that  the  intercept  is  ~  0.6  as  compared  to  an  in-line  orifice  (i.e. 
with  0  cross  velocity)  of  0.5  for  the  contraction  and  expansion  processes  only.  The 
difference  is  probably  related  to  the  friction  (-0.01)  and  turbulence  losses  which  at  this 
velocity  would  be  expected  to  be  small. 
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Figure  12  -  Correlation  of  Kl  for  Differing  Velocity  Ratio  and  Turning  Angle 


Table  4  -  Correlation  Coefficients  for  Non-Cavitation  Turbulent  Flow 


KL  A  j 

V  \  V  \ 

ft) B  ft)e 

ANGLE 

60 

75 

90 

105 

120 

A 

1.3353 

1.0708 

1.028 

0.8496 

0.7143 

B 

-1.1073 

-0.2061 

0.0638 

0.7475 

1.106 

C 

0.6787 

0.5241 

0.6118 

0.5415 

0.5592 

R2 

0.9915 

0.9877 

0.9899 

0.9950 

0.9956 

The  results  from  this  study  were  also  compared  with  pipe  bending  loss  coefficients 
typically  published  in  handbooks  and  fluid  dynamics  tests.  Typical  examples  are: 

45°  Kl=  0.42 

90°  Kl=  .090  rounded  bend  and  1 . 1  for  R/D  =  0 

180°  Kl=  2.2 

It  is  assumed  that  for  constant  area  bends  the  bending  loss  is  independent  of  Reynolds 
number.  In  fact,  no  discussion  is  found  to  suggest  that  KL  is  not  a  constant.  Since  for 
our  study  KL  varies  not  only  with  bending  angle  but  also  with  manifold  to  orifice  velocity 
ratio,  the  results  can  not  be  directly  compared.  The  comparison  that  was  made  was  to 
compare  the  value  of  Kl  for  various  V1/V2  ratios.  This  was  accomplished  by  first 
plotting  KL  vs  angle  at  constant  velocity  ratio  then  taking  the  average  slope.  The  results 
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are  plotted  in  Figure  13.  A  best  fit  linear  line  was  used  to  show  the  trend.  Note  that  as 
the  velocity  ratio  is  lowered  the  slope  and  value  of  Kl  approaches  that  of  the  constant 
area  bending  loss  coefficient.  It  would  be  expected  that  as  the  contraction  and  expansion 
losses  are  reduced  the  loss  coefficient  should  approach  that  of  a  constant  area  bending 
value. 


Turning  Angle  Impact  On  Loss 


Angle  ~  deg 


Figure  13  -  Comparison  of  KL  between  This  Study  and  Constant  Area  Bending 

The  application  of  our  results  to  predict  orifice  Cd  is  rather  straight  forward.  Obviously 
the  first  step  is  similar  to  that  described  above  utilizing  the  design  information  to 
determine  that  it  does  not  fall  within  the  cavitation  regime.  Once  this  is  determined  then 
Kl  is  determined  from  Figure  12  or  the  Equation  and  constants  provided  in  Table  4,  the 
head  loss,  HL,  is  determined  from  Equation  9.  Lastly,  Cd  is  then  calculated  from 
Equation  2. 

Compound  Angle  Direction  Feed 

Impact  of  Orifice  Variables  on  Kl  -  The  head  loss  coefficient,  KL,  for  the  compound 
angle  direction  feed  data  was  also  determined  as  a  function  the  manifold- to-orifice 
velocity  ratio.  The  results  are  plotted  in  Figure  14  and  compared  with  the  75  and  105 
degree  single  direction  results.  As  for  the  cavitation  data  the  turbulent  non-cavitation 
data  also  plot  on  a  single  curve  for  both  the  60_90  and  75_90  orifices.  However,  for  the 
non-cavitation  regime  the  comparison  suggests  that  while  the  90°  turning  angle  appears  to 
have  the  largest  impact  on  head  loss  the  orifice  turning  angle  impacts  can  not  be  ignored. 
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Figure  14  -  Velocity  Ratio  Impact  on  Kl  for  Compound  Angle  Orifices 


CONCLUSIONS 

The  results  from  this  study  demonstrate  that  the  simple  first  order  relationship  between 
Cd  and  Kcav1/2  as  well  as  HL  and  KL  are  valid  for  non-axial  orifices  with  cross  flow  and 
well  as  in-line  orifices  where  cross  velocity  is  0.  Also,  correlations  provided  should  aid 
designers  in  determining  where  cavitation  will  occur  and  define  the  Cc  in  addition  to  the 
head  loss  in  the  non-cavitation  regime.  Additional  test  efforts  should  be  directed  to 
extend  the  range  of  results  to  include  L/D  variation  for  all  angles,  inlet  r/R,  increased 
orifice  angles  for  compound  angles,  and  a  larger  range  in  orifice  diameters  as  well  as 
other  fluids  to  define  Reynolds  number  impacts.  It  is  expected  that  the  first  order 
relationship  will  not  be  adequate  for  larger  area  ratios. 

An  attempt  was  made  to  reduce  the  individual  equations  for  turning  angle  in  both  the 
cavitation  and  non-cavitation  regime  to  a  single  equation.  However  the  nonlinear  nature 
of  the  equations  did  not  produce  a  correlation  with  acceptable  error. 

The  significant  limitation  for  this  study  is  that  the  manifold  to  orifice  area  ratio  was  not 
varied  sufficiently  to  determine  its  impact.  Therefore,  the  resulting  correlations  are  valid 
only  where  the  area  ratio  is  sufficiently  small  that  the  term  can  be  ignored. 

The  conclusions  that  can  be  drawn  from  the  results  presented  are: 

CAVITATION  REGIME 

•  In  the  cavitation  regime,  Cc  is  controlled  by  the  cavitation  number 


18 


•  When  only  part  of  the  flow  enters  the  orifice  (i.e.  cross  velocity  flow)  the 
contraction  coefficient  is  governed  by  both  the  turning  angle  and  the  velocity  ratio 
V1/V2 

•  Inception  of  cavitation  occurs  at  a  Kcav  of  1.8 

•  Full  Cavitation  occurs  at  a  Kcav  of  ~  1.5  to  1.8 

•  For  the  compound  angle  orifices  Cc  is  more  impacted  by  the  initial  90°  turning 
angle  than  the  orifice  turning  angle 

NON-CAVITATION  REGIME 

•  In  the  non-cavitation  regime  Cd  is  likely  controlled  by  Reynolds  number  although 
in  this  study  Re  was  only  varied  by  velocity. 

•  The  head  loss  coefficient  Kl  when  only  part  of  the  manifold  flow  enters  the 
orifice  was  found  to  be  a  function  of  both  turning  angle  and  velocity  ratio  V 1  / Vc 

•  For  the  compound  angle  orifices  Kl  is  more  impacted  by  the  initial  90°  turning 
angle  than  the  orifice  turning  angle 


NOMENCLATURE 


SYMBOLS 

A  Area  in2 

Cc  Contraction  Coefficient,  A2/Ai 

Cd  Discharge  Coefficient 

Hl  Head  Loss,  lb/in2 

K  Loess  coefficient 

Kcav  Cavitation  Number 

P  Pressure,  psi,  psia 

Re  Reynolds  Number 

V  Velocity,  ft/sec 

W  Flow  Rate,  lb/sec 

GREEK 

p  Liquid  Density,  lb/ft3 

6  Turning  Angle,  degree 

SUBSCRIPTS 

1  Manifold 

2  Orifice  Exit 

a  Contraction 

c  Contraction 

out  Flow  rate,  velocity,  and  pressure  in  the  manifold 
turb  Turbulence 

f  Friction 
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